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We propose to describe the time evolution of quasi-stationary fluctuations near QCD critical point 
by a system of stochastic Boltzmann-Langevin-Vlasov-type equations. We derive the equations and 
study the system analytically in the linearized regime. Known results for equilibrium stationary 
fluctuations as well as the critical scaling of diffusion coefficient are reproduced. We apply the 
approach to the long-standing question of the fate of the critical point fluctuations during the 
hadronic rescattering stage of the heavy- ion collision after chemical freezeout. We find that if 
conserved particle number fluctuations survive the rescattering, so do, under a certain additional 
condition, the fluctuations of non-conserved quantities, such as mean transverse momentum. We 
derive a simple analytical formula for the magnitude of this "memory" effect. 



I. INTRODUCTION 

Mapping the QCD phase diagram as a function of 
temperature T and baryochemical potential fig is one 
of the fundamental goals of heavy-ion collision experi- 
ments. QCD critical point is a distinct singular feature 
of the phase diagram. It is a ubiquitous property of QCD 
models based on the chiral symmetry breaking dynam- 
ics 

0, 1 (see i for review and further references). Lo- 
cating the point using first-principle lattice calculations 
is a formidable challenge Recent progress and re- 

sults are encouraging, but much work needs to be done 
to understand and constrain systematic errors (see, e.g., 
Refs.[§-[ll| and reviews [llllllj for further references and 
discussion) . 

If the critical point is located in the region accessible to 
heavy-ion collision experiments it can be discovered ex- 
perimentally. The search for the critical point is planned 
at the Relativistic Heavy Ion Collider (RHIC) at BNL, 
the Super Proton Synchrotron (SPS) at CERN, the fu- 
ture Facility for Antiproton and Ion Research (FAIR) at 
GSI, and Nuclotron-based Ion Collider Facility (NICA) 
in Dubna 

The characteristic feature of a critical point is the 
increase and divergence of fluctuations. The non- 
monotonous behavior of event-by-event fluctuations, 
measured in heavy-ion collisions, as a function of the 
initial collision energy is a signature of the QCD criti- 
cal point [3 The estimates of the magnitude of 
the fluctuations in [lj| were based on the assumption of 
thermodynamic equilibrium, which is a reasonable first 
approximation at freezeout. For such stationary fluctu- 
ations the probability of a given value of a fluctuating 
variable is proportional to the exponential of the entropy, 
i.e., to the number of microscopic states with that value 
of the variable [13, [SI . 

In a dynamic environment of a heavy ion collision, 
the system continuously evolves with time. As long as 
the evolution is slow enough compared to the typical re- 
equilibration time, one can consider fluctuations as sim- 
ply tracking the evolving equilibrium conditions. How- 
ever, some fluctuating modes can be slower. In fact, it is 
precisely these slow modes which are of primary interest 



to us. These include fluctuations of conserved quantities 
and, most importantly, the critical fluctuations of the 
order parameter field a at the critical point. Fluctua- 
tions must keep readjusting to the continuously drifting 
equilibrium value. Can this quasi-stationary dynamics of 
fluctuations be described quantitatively? The purpose of 
this paper is to achieve this. 

We derive stochastic equations for the particle dis- 
tribution functions as well as the critical mode using 
fluctuation-dissipation relation in Section Mil We de- 
termine the corresponding equation for the correlators of 
the fluctuations in Section [V] and study its solution in 
Sections EII1 rVTHl andHXl Finally, in Section as an 
example of the application, we answer analytically, in an 
idealized regime, the long-standing question of the fate of 
fluctuations during the hadronic rescattering phase. We 
discover a "memory" effect, which protects not only fluc- 
tuations of conserved quantities. Notations introduced 
throughout the paper are indexed in Appendix [A] 



II. COMPARISON WITH RELATED WORK 

The time evolution of fluctuations has been consid- 
ered previously using different methods and/or in differ- 
ent contexts. Below we review some of this work in order 
to point out the new ingredients as well as the results of 
our approach. 

Quasi-stationary dynamics of fluctuations motivated 
Ref. [2^ ]. The relaxation of the correlation length, as a 
proxy to the magnitude of fluctuations, was studied using 
a model equation. Here we shall address evolution of 
fluctuations on a microscopic level and directly in terms 
of observable quantities. 

The quasi-stationary dynamics plays an essential role 
in the anomalous suppression of charge fluctuations, 
which has been proposed as a signature of the quark- 
gluon plasma formation at early times in Refs. [23|, l24j . 
and more quantitatively analyzed in [25[ (see also review 
Ref. [H[). The evolution follows a diffusion- type equa- 
tion, which means that fluctuations of larger spatial ex- 
tent relax slower. Therefore, as the size of the acceptance 
window is increased, the "memory" of the fluctuations 
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goes further back in time, allowing to probe earlier stages 
of the fireball evolution. 

The most interesting and not easily anticipated result 
of the approach we introduce here is the following. Al- 
though the "memory" effect is due to the slowness of the 
conserved charge fluctuations, the fluctuations of other 
quantities are also affected. For example, we show that 
the fluctuations of observables such as, e.g., mean trans- 
verse momentum px in the event also "remember" their 
earlier value. More precisely, if the chemical freezeout 
(the freezeout of inelastic reactions) has occurred near 
the critical point, the elastic collisions during the subse- 
quent evolution of the fireball do not completely "wash 
out" the critical point contribution to the mean p? fluc- 
tuations even on the time scales longer than typical colli- 
sional relaxation time. We can determine the magnitude 
of that effect by studying the microscopic nature and 
evolution of the slowest mode of fluctuations. 

The evolution of fluctuations in the vicinity of the crit- 
ical point has been studied numerically in Ref . [27j • The 
fluctuations were introduced by randomization of initial 
conditions, while the subsequent evolution was determin- 
istic. The essential ingredient of the approach we intro- 
duce here is the full treatment of fluctuations. I.e., fluc- 
tuations are driven by a random external source, acting 
at all times. The strength of the source is determined by 
fluctuation-dissipation relation. 

Stochastic Langevin-type equations were used to study 
hydrodynamic fluctuations near the critical point in [28j j . 
However, experimental observables such as, e.g., mean 
Pt fluctuations, are not directly related to hydrody- 
namic variables. The new ingredient in the present ap- 
proach is the use of kinetic Boltzmann equation. The 
degrees of freedom here are particle distribution func- 
tions, which directly translate into observable fluctuation 
measures. Strictly speaking, our approach lacks rigorous 
consistency of low-energy hydrodynamic description, and 
should be considered as only a model of the late hadronic 
stage of the heavy-ion collision. However, this relatively 
minor compromise allows us to address directly experi- 
mental fluctuation measures and make quantitative pre- 
dictions, rather than limiting the study to density fluc- 
tuations. As a test of the new approach we shall derive 
some results of Ref. [H| in Section IVHTBl 

In a different context, the stochastic Boltzmann- Vlasov 
type equations have been used to estimate the rate for 
hot electroweak baryon number violation [29l - l3~l| . 



III. THE FORMALISM 

To model the fireball evolving through the phase dia- 
gram near the critical point, we consider a relativistic sys- 
tem of particles interacting with a scalar field <r, coupled 
to a thermal bath at temperature T. We are considering 
late hadronic phase of the fireball expansion. The par- 
ticle density is assumed to be already sufficiently small, 
so that their motion can be considered classically, using 



Boltzmann equation. The mass m of the scalar field, 
on the other hand, is considered to be sufficiently small 
compared to 2irT, so that the field can be treated clas- 
sically, using field equations. This condition is fulfilled 
sufficiently close to the critical point. 

A. Equations of motion 

Since the field a is a Lorentz scalar, we assume that 
the coupling of it to the particles affects their mass (as 
opposed to, e.g., chemical potential, which is a Lorentz 
vector). As an example, one can consider coupling of 
nucleons aNN or pions aim to the a field in the chiral 
sigma model. Thus we are led to consider classical motion 
of particles with variable mass M(a) [32J, which depends 
on the local value of the scalar field a. The action of the 
system is given by: 

S = fd 3 x^(d fJl ad^a-U{a)) - JdsM(a), (1) 

where the last integral is taken over the worldline of a 
particle with variable mass. The corresponding equations 
of motion are given by: 

d 2 a + dU/da + J ds (dM/da) = 0; (2) 

dp>*/dT = d»M(a), with p^ = Mdx t "/dr, (3) 

where M = M(a) is the local value of the variable par- 
ticle mass. One can check that the motion governed by 
equations ([3]) preserves p^p^ — M(a) 2 = along the par- 
ticle trajectory. 

The Boltzmann equation for the distribution function 
f(x,p) of such particles in the external field a reads [33j j : 

or in a more physically transparent, non-covariant, form: 

f + v-Vf- (VM/ 7 ) • (df/dp) + C[f\h = 0, (5) 
where 

v = p/i^M) and 7 = (1 - v 2 )~ 1/2 (6) 

is the particle velocity and relativistic gamma-factor re- 
spectively. 

The collision integral C [/] in Eq. (j4j gives the collision 
frequency for all particles with momentum p (near space- 
time point x) in the rest frame of those particles, while 
C /"f is that frequency in the lab frame. Using equations 
of motion @ one can show that the Boltzmann equation 
(HI) (or |5])). implies continuity equation for the particle 
number current 

d,r + I C[f\h = , where f = [ /p"/(M 7 ), 

(7) 



3 



with 



where 



(2tt) 3 



(8) 



I.e., particle number can only be changed (if at all) by 
collisions. 

The equation of motion for the scalar field a in the 
presence of particles with distribution /(ir,p), following 
Eq. ^ , is given by 



d 2 a + dU/da+(dM/da) / f/j = 0. 



(9) 



Coupled equations ((U) and © describe evolution of the 
particle distribution / and the scalar field a. These equa- 
tions are conceptually similar to Vlasov equations in elec- 
trodynamics. The difference is that the classical field a is 
a Lorentz scalar. There is also certain limited similarity 
with the nuclear mean-field approach 34j . 

Our goal is to extend the above formalism to the de- 
scription of fluctuations in the system. In application 
to linearized Boltzmann equation this has been done by 
Fox and Uhlcnbeck and others f35U38j . Here we shall ex- 
tend this formalism to linearized Boltzmann- Vlasov type 
coupled equations (0| and ©■ 



B. Linearized equations 

For a given constant field a, the Boltzmann equation 
has a stationary solution, which is also constant in space, 
f a {p), satisfying C[f a ] = 0. This is Boltzmann distribu- 
tion for particles of mass M(a) at arbitrary values of 
temperature T and chemical potential (i: 



U{p) = e 



= pM/T e -7(p)M/T _ 



(10) 



The values of T and /i depend on the total particle num- 
ber (if it is conserved by collisions) and on total energy, 
if the system is closed, or by conditions of equilibrium 
with the thermal bath if it is open. 

The equilibrium value of a is determined by 



dU/da + (dM/da) / f a /j = 



(11) 



where the second term can be viewed diagrammatically 
as the contribution of a thermal tadpole. 

We linearize the equations for a and / by expanding 
around their equilibrium value. The deviation of / from 
its equilibrium value fa- will be parametrized, as usual, 
by function h: 

f = f*0- + h). (12) 
The linearized Boltzmann equation then reads 

h-ag/(jT)+v Vh+l[h] = 0, (13) 



g = dM/da (14) 

and I[h] is the linearized collision integral: 

C[f] = 1 f rJ l[h]+0(h 2 ). (15) 

Note, that both C and I depend on the local value of the 
field a (through the dependence of the particle mass M), 
and we used the property of the equilibrium distribution 

c[/ CT ]=o. ~ ' ; 

Shifting the notation for a so that a = is the equi- 
librium value (solution of Eq. (JXTJ) ) we can write the lin- 
earized equation @ as 



a - V 2 cr + m 2 a + g f h/j = . 



(16) 



where we defined the "in-medium" mass m of the the a 
field quanta as 



d 

do 



(17) 



cr=0 



with ml = d 2 U(0)/da 2 . The last term in Eq. (TTT]) can 
be recognized as the one-loop thermal contribution to the 
vacuum mass m . 

In the system we considered so far the dissipation (en- 
tropy increase) is entirely due to the collision term C [/] . 
In a more general, and more realistic, case when the field 
a interacts with other particles in a heat bath, one can 
describe the additional dissipation effects adding a term 
T <T into the l.h.s. of Eq. (ITB1) . 



C. Noise and its correlators 

Equations (O and ([9]) describe evolution of functions 
/ and a averaged over the time scale of many particle 
collisions. Fluctuations of / and a can be characterized 
by correlation functions (also averaged over many colli- 
sions). 

In order to describe these fluctuations we follow the ap- 
proach of Ref. (36| and introduce random noise terms. We 
shall determine the correlation functions of these noise 
terms following Refs. [H, HI, HI, 53] in the linearized 
regime. For that purpose we shall cast equations in the 
following first-order form: 



h - irg/(-fT) 
n + r 7r — 
a — 7T = 0. 



v- Vh+X[h] = £; 



vV 



m 2 a 



9 I /oV7 



m 



(18a) 
(18b) 
(18c) 



where we introduced noises £ and rj and, for generality, 
the additional dissipation term Tqtt. It is worth pointing 
out that the field a is stochastic with our without the 
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noise rj because of the couping to the particles (last term 
on the l.h.s. of Eq. (pEt ). 

To determine the correlators of the noises £ and r/, 
we use fluctuation-dissipation relation. The probability 
distribution of the fluctuating degrees of freedom [2CJ, [21| 



V[f,a]~expS[f,a] 
is determined by the entropy function 



(19) 



S[f,<r] 



d 6 x 



t It 



(Va) 2 



/(log/-l) 



+ [/(a) + / (M(a) 7 - M )/ 



(20) 



where T is the temperature of the external heat bath. 
The first term is the well-known Boltzmann entropy (H- 
function), while the second term is simply — {E — fiN)/T, 
where E is the energy of the system and fi is the chem- 
ical potential. This term is the contribution of external 
reservoir to the (fluctuations of) entropy. 

In the linear approximation, we can consider noise 
to be Gaussian, and all non-trivial information about 
it to be in the correlators such as (£(xi,pi)t;(x 2 ,p 2 )}, 
(r](xi)r](x 2 )) and (£(x\,p)r](x 2 )}. To determine these 
correlators we expand the entropy to quadratic order 
(note cancellation of terms linear in h due to log f a — 
(tx-M(a) 7 )/T) 



S< 2 



where 



d s x 



foh 2 



1 

T 



2„2\ 



(Ver) 2 + m 2 a 



m 2 = mi-T 



d 2 U 
da 2 



(21) 



(22) 



is the same "in-medium" mass of the a field quanta al- 
ready defined in Eq. (|17[) . as can be verified by using 
Eq. ([ID]). 

We follow refs. [H, HI, HH to define the "entropy ma- 
trix" (or, more precisely, operator) E: 









-I 







foh 
ir/T 

(-V 2 a + m 2 a) /T j 



so that Eq. (j2"Tj) can be written as 

S {2) = -h-Eh, 
2 



where h = (h, n, er) denotes the (infinitely dimensional) 
"vector" whose components are the degrees of freedom 
of the system, with the scalar product defined as 



(h,7r,a)-(h',7r',a') 



x 



hti 



(25) 



Similarly, equations (|18l) can be also cast in matrix 
(operator) form 



where 









-1 







/» ::/) c 



--Kg/{iT) + v Vh + l[h] 



(26) 



(27) 



Then the correlator of the noises, combined into a vec- 
tor £ = (£, 7/,0), can be expressed in terms of the ma- 
trix/operator Q defined as 



(€(*l)®€(*2)) =2QS(t 1 -t 2 ) 



(28) 



and given by the usual fluctuation-dissipation relation 
(see Refs. jaESHi) 

2Q = GE" 1 +E~ 1 G t . (29) 

Using Eqs. (p3]l. ([2T]l and ([29j one can now find Q: 









-( 







I[h/f ]+tf[h]/tf 
2T q Ttt 




where we defined operator K. as 

I[h] =K[foh] 



'(/C + £+)[/i] N 
2T Tir 




(30) 



(31) 



One can show [40J that the operator K is self-adjoint for 
elastic collisions, but we leave equations in a more general 
form. Equation (|30|) together with Eq. ([28)1 translates 
into the following explicit expression for the correlators: 



= {JC+ /C t )(2^) 3 ,5 3 (pi - P2 )5\x 1 - x 2 ); (32a) 



{■q{xi)r](x 2 )) = 2r T<5 4 (a;i - x 2 ) 
(^{xi,Pi)v(x 2 )) = 0. 



(32b) 
(32c) 



It is easy to reco gniz e in eq. (|32a[) the generalization of 



the result of Rcf. 



, . One can also observe that the in- 

(23) teraction between the particles and the field a does not 
manifest itself in any modification of the corresponding 
noises. That should be expected given the physical origin 
of the noise: collisions and the interaction with the ex- 
ternal reservoir. The correlations are local in coordinate 

(24) space, which also correctly reflects their origin. 



IV. STATIONARY, EQUILIBRIUM 
FLUCTUATIONS 



The quantity directly accessible by experimental mea- 
surement is a two-particle correlator (Sf(r\Sf( 2 \). Before 
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we begin studying time evolution of fluctuations let us de- 
rive the stationary, equilibrium value of the two-particle 
correlator and compare with existing results. 
To linear order in fluctuations, 

/ = f a (l + h) = / (1 +h- gal OT)) + 0(a 2 ). (33) 

The equal-time correlators of h and a are contained in 
the matrix elements of the correlator of h which in equi- 
librium are given by 



(h®h) = E~ 1 



(34) 



according to 1)240. Using the explicit expression (|23)) for 
the components of the entropy matrix we thus find 

(Sf{i)Sf(2)) = (/(i)/(2)) - /o(i)/o(2) 

= /0(1)/0(2) ((^^) (1) { k -£) W ) 



— /o(l)<5(l,2) + 



g 2 /0(l)/0(2) 



7(1)7(2) 



T 



(35) 



where subscripts (1) and (2) refer to the points in the 
phase space (xi,pi) and (x2,P2) where the scripted 



quantities are to be evaluated, Sn,2) 



S 3 ( Xl 



x 2 )(2tt) 3 5 3 { Pi ~ P2 ) and D {1 , 2) = (- V 2 + m 2 )- 1 S 3 (x 1 - 
aj 2 ). Integrating over i 12 one obtains the known result 
for the equilibrium fluctuations of particles coupled to 
classical scalar field a pjl, |4l| : 



V 1 {8vp&v p >) = fo6pp> + 



fo fo 



m 2 T 7 7' 



(36) 



where V = j d 3 x is the volume, S pp r = (2ir) 3 S 3 (p — p'), 
and we denoted the momentum space distribution as 

v p = J d 3 x f(x,p), such that J v p = N (37) 
is the total number of the particles. 



V. TIME EVOLUTION OF FLUCTUATIONS 

Now we want to consider the time evolution of fluc- 
tuations as the parameters of the system, most impor- 
tantly m, change. 

If we take the initial probability distribution for fluctu- 
ating variables to be Gaussian, in a linear system the fluc- 
tuations will remain Gaussian at all times. This can be 
verified directly, by converting the generalized Langevin 
equations (|26[) into corresponding Fokker-Plank equation 
for the probability distribution P[h; t]: 



d ( dP 
P =dH-( GhP + ®3h 



(38) 



Parametrizing the probability using (time-dependent) 
operator S: 



P = (detEr 1/2 exp 



1 



(39) 



and substituting into Eq. (|38[) . one finds equation for S: 
± = -GS - £G f + 2Q. (40) 

Note that S = E _1 is a stationary solution of this equa- 
tion, as expected from Eq. (|34"|) and the fact that 

(41) 



(h®h) = S. 



This verifies the fluctuation-dissipation relation (1291) . 

The equation (|4"0|) for (d/dt)(h <£> h) can be also de- 
rived directly, by applying Eqs. (|2^|) and (|28p . without 
assuming Gaussianity (|39[) . 

Equation (|40p can be formally integrated from initial 
time, taken to be 0, to arbitrary time t: 

E(t) = V(t,0)S(0)V t (t,0) + 2 f dt'Y(t,t')QV^(t,t'), 

Jo 

(42) 

where evolution operator V(t,t') satisfies 



V = -G(t)V, 



V(*,0=Te«p (-J 



dt" G(0 



(43) 



(44) 



Equation (|40|) can be now used to study the time evolu- 
tion of the fluctuations, provided, of course, they remain 
small, so that linear approximation is valid at all relevant 
times. 

A more useful equivalent form of equation (|40[) is 



^ = -GE-EGt-« 
dt dt 



where S 2J ... 1 . 

.(45) 

In this form, and for E = const, it describes relaxation of 
X to its equilibrium value E _1 . The solution is given by 

E(t) - V(i, 0)E(0)Vt(i, 0)- J dt' V(t, O^^^Ct, *'). 

(46) 

which is equivalent to Eq. (l4"2"l) . 

Determining the time-dependence more explicitly in 
the general case is a complicated task, and should per- 
haps be part of numerical modeling of a more realistic 
system. As an illustration of the use of Eq. (l4"6l we shall 
address the important question of the fate of the critical 
point fluctuations after chemical freezeout. To prepare 
for this, we shall briefly discuss the way conservation of 
particle number is reflected in the equations (Section lVIl) . 
and then analyze the evolution of fluctuations for time- 
independent G and E (Sections I VIII and IVIII|) . 
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VI. CONSERVATION OF PARTICLE NUMBER 
AND FLUCTUATIONS 

Chemical freezeout is a moment in the history of a 
heavy-ion collision fireball, when inelastic reactions be- 
come too infrequent to modify the chemical composition 
of the system. In other words, the number of particles of 
a given species is conserved during subsequent evolution. 

To model the evolution past the chemical freezeout, we 
shall require that the collision integral C [/] conserves the 
particle number. According to Eq. ([7]) this requires 

/c[/]/ 7 = (47) 
J P 

to be valid for all /. The linearized collision operator I 
in Eq. (Tl"5l) therefore obeys, for all h, 

M[l[h}}=0 (48) 

For future convenience, we have introduced notation 

M[h] = ^ (49) 
J P Jo 

for the average over equilibrium distribution /o- Impos- 
ing the condition that operator JC defined in Eq. (J2J) is 
self-adjoint [40(, one can see that Eq. (|4"5|) implies 

Z[const] = 0. (50) 

In other words, operator X has a zero mode. This is 
also evident from the fact that constant h corresponds 
to changing the value of /j, in the equilibrium distribu- 
tion (fT0|) . and that C[f a ] = for arbitrary /z. 



VII. SOLVING EVOLUTION EQUATIONS 
WITH CONSTANT COEFFICIENTS 

The evolution operator V in Eq. (|44l) for the system 
with time-independent G can be written in the form: 

Y(t,t')=J2^ Mt ~ t ' ) hx^e x , (51) 

A 

where the sum goes over all solutions of the following 
eigenvalue system 

Xh x = Gh x , (52) 

and vectors 6\ form the dual (adjoint) basis with respect 
to the one formed by vectors h\, i.e., 

Ox - /'a <V (53) 

If the system is also spatially homogeneous, it is conve- 
nient to apply Fourier transformation with respect to the 



space coordinate x to equations (|52|) . Given the defini- 
tion of operator G in Eq. pT)) . we find: 

-\h-ng/(~/T) + ivqh + l[h] = 0; (54a) 
-Att + r 7T + q 2 a + m 2 a + gn M[h/j] = 0; (54b) 
-Act - tt = 0; (54c) 

where we defined 

«o = / fo (55) 
J P 

- the equilibrium density of the particles. 

VIII. THE SLOWEST MODE 

The slowest mode corresponding to the lowest eigen- 
value of the eigensystem (|5~4f is of primary interest to us. 
In this section we shall determine it. 



A. Zero mode 

We begin by considering the simpler case q = 0. Since 
operator I has a zero mode h = const (|50|) , let us separate 
it by writing 

h = h + h, where h = M[h]. (56) 

The zero eigenvalue Ao = of Eqs. ([54|) corresponds 
to the solution such that h = n = (so that Eq. (|54a|) is 
trivial) and 

h Xo : gn M[lh]h = -m 2 a (A = 0). (57) 

Note that this mode is predominantly a near the critical 
point (i.e, h — > as m — > 0). This may appear surprising, 
since the mode Ao is due to the particle number conser- 
vation, while a is not a density of a conserved quantity. 
However, near the critical point the fluctuations of par- 
ticle number density are dominated by their mixing with 
a [Elm, whose fluctuations diverge. Eq. (f57l> also shows 
that in the limit g — > at fixed m, the Ao mode is pre- 
dominantly h, as it should be if a is decoupled. 

In order to find the dual vector 9\ we need some in- 
formation about all other modes h\, since B\ must be 
orthogonal to them (|53"|) . Applying operator M. to equa- 
tion (|54a[) and using (|48|) . we find, at q = 0: 

h x : h = ^M[lh\a (A^0). (58) 

This equation contains all the information about the non- 
zero modes that we need to determine 9\ . 

In order to simplify subsequent linear algebra manip- 
ulations, we shall define two convenient basis vectors: 

h: h = ir = a = and h = gM[l/~f]/T; (59) 
<t: h = 7r = h = and a = 1. (60) 
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In terms of these vectors, Eq. (|58l) means that, for any 
nonzero eigenvalue A ^ 0, eigenmodes are given by 

h\ = a + h + (terms with h = a = 0) , (61) 

while Eq. (l57l) for the zero mode can be written as 

h Xo = A 2 <r - m 2 h , (62) 

where we defined 



(63) 



Equations (jBTj) and (|62|) together determine the orienta- 
tion and the length of the dual vector 



0a„ = - 



1 <T-0 



V A 2 + m 2 ' 



(64) 



which satisfies the defining orthonormality condi- 
tions (f5"3")) . We defined another convenient vector, related 
to ET 



(65) 



such that • h = 1. One can also check that 9\ is the 
A = eigenvector of X0\ = &0\, as it should be. 



B. Hydrodynamic mode and diffusion coefficient 

The mode Ao, corresponding to conservation of the 
particle number, is hydrodynamic in the sense that, for 
small q, An — 0(q 2 ). The ratio Xo/q 2 — D defines the 
corresponding diffusion coefficient D, which we can ex- 
tract from equations (f54|) . 

As we did deriving Eq. (1551) for nonzero modes at q = 
0, let us apply operator M. to Eq. (|54a[) . Now, at q ^ 0, 
we find 



Xh -iM[v- qh] = X±M [l/i\a . 



(66) 



We need now to express h in terms of h and a using 
equation (|54al) . Since h — » as q — > 0, one can see that 
h must begin at order q. Keeping in Eq. (|54a[) only terms 
of 0(q) we obtain: 



iv • qh+X[h] =0. 



(67) 



We would need to invert operator I to express h in terms 
of h. For generic operator I we shall define function 
ip(v 2 ), which solves the equation 



X\v ib] 



(68) 



The fact that solution can be found in this form follows 
from isotropy of the collision operator and equilibrium 



distribution function. In terms of ip, the solution to equa- 
tion ()67j) is given by 



-iq • v iph . 



(69) 



Substituting this into Eq. (|66| we find, instead of 
Eq. (J5J), 



h(\ -D q 2 ) = \o^M[l/~f}<J, 



T 



where we denoted by D 



D a = -M[v 2 ^] 



(70) 



(71) 



the diffusion coefficient for the particle gas with fixed 
mass (the limit g — > 0). Putting together Eq. ([70)1 
and (|57|) . which remains valid to the order in q 2 we need, 
we find 



An 



Dq 2 +0(q i ), where D 



A 2 



rA). (72) 



The fact that D — > as m 2 — > is to be expected on 
gener al g rounds from the hydrodynamic relation D = 
&/x [431 j where a is the conductivity, and x is t ne 
susceptibility of the particle number, and the fact that 
X ~ 1 /m 2 (2a , |43| . Within our microscopic approach: 



(SN 2 ) 
VT 



no A 2 
T 



according to Eq. (|93"1) which we encounter later. 



(73) 



IX. FASTER MODES 

This section is a slight detour from the main thread of 
the paper. We have already accumulated all information 
about the zero and even non-zero modes, Eqs. ([ST|) . 
that we need to study the "memory" effect in fluctua- 
tions (Section [Xj. However, it might still be interesting 
to look at the structure of the non-zero modes in more 
detail, to understand better the properties of the system 
of equations we are solving. 

So far we have not used any information about the lin- 
ear collision operator I beyond the conservation of the 
particle number and isotropy. For the sake of analytic 
transparency, and within this Section only, we shall as- 
sume here that all eigenvalues, but one, of the operator X 
are equal to the same value r , which has the meaning 
of an average relaxation rate. The exception is the zero 
eigenvalue, corresponding to the condition This ap- 
proximation is well known and is due to Refs. [44l. |45| (see 
also Ref. (H). 

Operator I should also respect the condition that op- 
erator JC, defined by (|3"Tj) is self-adjoint. All the above 
conditions are satisfied by the operator 



X[h] =T~ 1 {h-M[h}) 



(74) 
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Another simplification we shall adopt in this section 
only corresponds to assuming t _1 3> r ,m, i.e., that 
t -1 is much faster than any other rate in the problem. 

We emphasize that, although these approximations are 
physically sensible, they are only used here to make a 
transparent analytic treatment possible, illustrating the 
properties of the system we study. 

Substituting Eq. (|58|) back into (|54a|) and solving for h 
we find, at q — > 0, 

h= T ^§;{M[l/i\-l/<y)<T. (75) 

Now substituting h given by Eqs. ([56| . (|58]l and (JT5J) into 
Eq. (I54bp , we find the equation determining the eigenval- 
ues A: 

A 2 - Ar + m 2 + - 0. (76) 

1 — At 

where we used (|63l) and defined also 

Al^ 9 -^M[lh\ (77) 

Equation (f76|) has three roots. For the scale hierarchy 
we consider, t _1 3> To, to, there are two roots of order Tq 
or m and one root of order r . The two smaller roots, 
to leading order in r, satisfy the quadratic equation 

A 2 - Ar + to 2 = , (78) 

where 

to 2 = to 2 + A 2 ; (79) 

r = T + A 2 r; (80) 

and thus 

A lj2 = -r/2±iv/™ 2 -(r/2) 2 . (8i) 

At this point one can see that rh is the rest mass (pole 
mass) of the quasiparticle a. It is different from the 
(static) scree ning mass to, and does not vanish at the 
critical point flJIIS III , where m -> 0. 

Also, Eq. (|80|) shows that the full dissipation rate T 
contains contribution ApT from the interaction of a with 
the particles. In principle, one could start with To = 
(closed system) and consider the particle collisions to be 
the only source of the dissipation. 

The third eigenvalue is given, to the leading nontrivial 
order in r, by 

A 3 = r- 1 + (A 2 - A 2 )r. (82) 

(Applying Cauchy-Bunyakovsky-Schwarz inequality to 
Eqs. ([63]) and (J77J), one can see that A < A r .) 
Finally, any function h± satisfying 

M[h±]=0 and M[h±/i\ = Q (83) 



solves the eigensystem with a = ir = and A = r _1 . The 
linear (eigen) space defined by Eqs. (JH3J) is infinitely di- 
mensional and, correspondingly, A = r _1 is an infinitely 
degenerate eigenvalue. This degeneracy is not lifted be- 
cause modes h± do not mix with the modes correspond- 
ing to eigenvalues A^, i = 0,1,2,3. This is a convenient 
feature of the Anderson- Witting approximation (|74p . 

Returning to the Ao mode, with operator I in the form 
given by Eq. ([74]) . equation can be solved: — t, 
and D =M[v 2 }t/3. 



X. FROM CHEMICAL TO KINETIC 
FREEZEOUT 

A. Preliminaries 

As an application of the formalism, let us consider 
the following long-standing problem. Assuming that 
the chemical freezeout occurred near the critical point, 
how much of the fluctuation signal survives until kinetic 
freezeout? A more precise and detailed answer to this 
problem will likely require a numerical simulation. Here 
we want to illustrate the mechanism, and make a simple 
estimate of the effect. To that end we shall make several 
simplifying assumptions, in order to maintain analytical 
control. In essence, we shall assume that separation of 
different relaxation time scales is sufficiently large for us 
to be able to focus on only the most relevant modes. 

Chemical freezeout is characterized by "freezing" of 
inelastic reactions. This means that the number of each 
individual particles is conserved (particles in the same 
isospin multiplet could be considered as different inter- 
nal states of the same particle, to allow for quasielastic 
collisions). This, in turn, means that any measure of 
fluctuations of a conserved number of particles should 
not change. More precisely, it can only change by diffu- 
sion, which we shall assume here to be the slowest scale 
in the problem (i.e., we work in the q — > limit). 

On the other hand, fluctuations such as those of 
mean px , which is not a conserved quantity, must evolve 
between chemical and kinetic freezeout, at which point 
they are "frozen" and eventually observable. The form 
and the amount of this evolution we shall now discuss. 

During the interval between chemical and kinetic 
freezeout the typical time scale, r e , of the evolution of 
the system is much slower than the inverse elastic colli- 
sion rate r and the scales T _1 associated with the re- 
laxation of the a field. In a realistic heavy-ion colli- 
sion T e = O(10 — 20) fm (order of fireball size), while 
r _1 , t = O(0.5 — 2) fm (typical hadronic scales). Thus 
we shall assume r e 3> t, 

In order to be able to obtain analytically tangible solu- 
tion we shall take into account the effect of the change of 
only one parameter: to - the screening mass of a. Since 
the fluctuations are singular as 1/m 2 near the critical 
point, the effect of change of to could be assumed to be 
dominant, compared with the change of, e.g., equilibrium 
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distribution functions fo (e.g., via change of T), which we 
shall consider fixed, for simplicity. As a concrete exam- 
ple, one could consider evolution of m determined by the 
model in Ref. [22j]. As we shall see, the actual time de- 
pendence of m will not matter, as long as it is slow, which 
is helped by critical slowing down [221 ] . 

The physically reasonable assumptions spelled out 
above are needed to make the analytic results attainable 
and usefully transparent. These assumptions can be re- 
laxed, e.g., via a numerical simulation, at the expense of 
analytic control. Our main purpose here is to illustrate 
the mechanism in the most transparent way possible. 



B. Evolution of fluctuations 



prefactor is close to unity near the limit we have been 
working in: q — > (the fact that D ~ m 2 — > near 
the critical point also helps). In general, the importance 
of the prefactor depends on the size of the region over 
which the fluctuations are measured. Since 1/Ao is an 
estimate of the time, td, it takes for a fluctuation to dif- 
fuse over this region, the factor can be estimated roughly 
as exp[— (t — t')/TD\- Below we shall consider the case 
when the region is large enough, so that td 3> r e . 

We are now ready to apply Eq. (|4"6"]) . We shall take the 
initial time t — to be the time of chemical freezeout, 
and the final time t = t k the time of kinetic freezeout. 
At chemical freezeout, t — 0, the fluctuations are equi- 
librated and S = 0. Thus at kinetic freezeout, t = t k , 
equation (|4l)]) gives, upon integration, 



We begin by determining the evolution operator 
Y(t,t'). If we choose the interval, t — t' , so small that 
we could neglect the change of G (i.e., change of m) and 
consider it constant, then we could integrate Eq. (j44j) and 
obtain V given by Eq. ([ST]) . 

If the interval t — t' is also large compared to relaxation 
scales A" 1 , for all i ^ 0, only the term corresponding to 
the zero mode A — Ao — 0(q 2 ) will survive in Eq. ([ST]) : 

Y{t,t') = e" Ao( *~ t ' ) /i Ao (g)0 Ao -|-(exp. small terms) . (84) 

In order to extend this result to longer time intervals 
over which the change of G cannot be neglected, we use 
the property 



Y(t,t')=Y{t,t n )...Y{t 2 ,t 1 )Y{t 1 ,t') 



(85) 



and subdivide t — t' into smaller intervals satisfying 
A^ 1 <C t n — t n -i <C T e . Our assumption of scale hier- 
archy is needed to make such a choice possible. Using 
Eq. (gH) we then find 

V(M') =e-ti x ° dt h Xo (t)®8 Xo (t') 

x (0 Ao (t) • h Xo (t n )) . . . (e Xo (h) ■ h Xo (t')) . (86) 



In order to evaluate dot products in Eq. (1861) . we use 
explicit form of eigenvectors h\ and Q\ , given by 
Eqs. (J62J) , (|64]), and find, e.g., 



e Xo (h) ■ h Xo (t') = 



A 2 



m 2 (t'^ 



A 2 + m 2 (t 1 ) 



(87) 



There is a string of such factors in Eq. ([86]) and, multi- 
plying them successively, one finds that all but the first 
and the last factor A 2 + m 2 (t) cancel, leaving 



Y(t, t') = e- J? ^ fc Ao (t) e Xa [t') %\ m JQ 

= e -&\ dt h Xo (t)®e Xo (t). 



The only dependence on the initial time t' remains in the 
exponentially decaying prefactor. Since Aq = Dq 2 , this 



E(tk) = h Xo (t k ) ® Xo (t k ) 

x (E- 1 (< c )-E- 1 (i k )) Xo {t k )®h Xo (t k ). (89) 



A shorter way to derive Eq. (|89|) is to observe that the 
actual time-dependence of m(t) is not important, as long 
as it is faster than the diffusion: tjj r e (but r e ^ A^ 1 ). 
Choosing m(t) to have an (almost) instantaneous step 
from m c to m k , and constant at all other times one can 
then find solution (|89l) using equation (|46l) with initial 
condition S(0) = E _1 (t c ) - E" 1 ^). 

Taking into account Eq. (|23[) , which for convenience we 
write, using notations 1^ for the unit operator th[h] = h 
and 7T for the basis vector (h = 0, n = 1, a = 0), 



E-\t) 



/o~% 



Ttt (£) 7T - 



m 2 {t) 



(90) 



together with Eqs. (J62J and ([64]), we find for S = E" 1 - 
at kinetic freezeout time 



E(tk) =E" 1 (t k ) 



T 



T \ A 2 <r 



A 2 & 



m^h 



A 2 



A 2 



(91) 



where m Cj k = m(i c ,k) is the value of the a screening mass 
at chemical/kinetic freezeout. The last term, containing 
1/m 2 , is the "memory" effect, due to the freezing out of 
conserved particle number fluctuations. 



C. Two-particle correlator and "memory" 

In order to translate Eq. (|9"T]) into observed fluctua- 
tions, we should recall that S = (h <g> h) and apply 
Eqs. ([55]) (all but the last equality) to calculate the 2- 
particle correlator. The fluctuations of the momentum 
space distribution of particles (|37[) at kinetic freezeout 
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are thus given by 



D. Estimating the "memory" factor 



V 1 {5v p 8v p >) = foS pp i + 



fo fo 



miT 7 7' 



91 
T 



1 



/o/o 



, A 2 / 7 + m|M[l/ 7 ] A 2 /y + m^[l/ 7 ] 



A 2 



which should be compared to Eq. (|36[) with m = m^. To 
check that the effect of the additional term is to preserve 
the particle number (multiplicity) fluctuation at the value 
it attained at chemical freezeout, let us calculate that 
fluctuation by integrating Eq. (|92p over momenta p and 
p' and using Eq. (f3"T)l . Normalizing by the total number 
(N) = HqV for convenience, we find 



(N) 



= 1 



A 2 

2 

mi 



1 



1 



A 2 = 1 



A 2 

2 



in 



(93) 



where we used definition ([6"3"]1. We see that, as expected, 
the effect of the "memory" term is to keep multiplicity 
fluctuations from changing after chemical freezeout. 

The effect which is less obvious is that the "memory" 
term also contributes to fluctuations of quantities which 
are not conserved. We shall keep discussion as general 
as possible, but to be less abstract, we shall consider 
fluctuations of mean transverse momentum px per event, 
which is one of the most common "intensive" measures 
of fluctuations. This fluctuation can be also expressed 
via the correlator ([92} (see, e.g., Ref. [HI): 



((Spt) 2 ) 



1 



(AO 2 V P . 

where we defined 



(pr -Pt)(p' t ~ Pt){5v p 5v p >) , 

(94) 



p T = M [p T ] ■ 



(95) 



Normalizing by (N) to remove trivial system-size scaling, 
we find 



(N)((5 Pt ) 2 )=M[(pt-Pt) 2 ] 
+ i^M[{p T -PT)hf 



1 - r m 



where we introduced 



(96) 



(97) 



Eq. (|96p shows that the critical contribution 0(l/m 2 ) 
can, under certain conditions, survive through the 
hadronic rescattering stage until kinetic freezeout. Com- 
pared to the value at chemical freezeout, the C(l/m 2 ) 
term is attenuated by the factor r m (l9"7T) which, if the a 
screening mass at kinetic freezeout, m^, is of order A or 
smaller, is a non-negligible fraction of unity. 



Let us now estimate the "memory" factor (|9"T|) . The 
value of r m depends quite strongly on the ratio of to 
A. For fluctuations to survive, m^/ A cannot be large. 

The estimate for A can be made using Eq. (|531 . In or- 
der to do this correctly we need to generalize our analysis 
to include more than one species of particles: nucleons (2 
spin and 2 isospin states), pions, etc. We then find that 
the expression for r m in equation (|97[) still holds, with 
A 2 receiving contributions from all species: 



A = A 



nucleons 



A ■ 

p IO IIS 



(98) 



Choosing, for example, top SPS energy freezeout con- 
ditions T = 168 MeV and [i B = 266 MeV [H, 
we find for the contribution of nucleons A nuc i cons ps 
430. (ffp/10.) MeV. We take g p w m p J U ~ 10 as an 
estimate of the coupling of a to protons. 

The estimate for the contribution of pions is A p i ons w 
110. (g n /2.) MeV, where g n s» Gjm^ ~ 2., using the esti- 
mate for G from Ref. [l!| . The estimates for the contribu- 
tion of antinucleons and kaons are similarly small, com- 
pared to A nuc ieons- Summation in quadratures increases 
the estimate for A by less than 10% over A nuc i oons : 
A « 460 MeV. 

Thus, at top SPS energy, the critical pr fluctuations 
survive at least half as well as the particle multiplicity 
fluctuations (r m > 1/2) until kinetic freezeout, if the a 
screening mass at the freezeout does not exceed mt < 
"Ta « 300 MeV. 



XI. SUMMARY, DISCUSSION AND OUTLOOK 

In summary, we introduced an approach to studying 
time-dependent quasi-stationary fluctuations near QCD 
critical point by combining stochastic Boltzmann equa- 
tion with an equation of motion for a scalar field, describ- 
ing the "soft" critical mode. We obtained the general so- 
lution of the linearized system and studied its relaxation 
modes. We focused on the slowest (diffusion) mode and 
analyzed its effect on the evolution of fluctuations after 
chemical freezeout. 

One of the consequences of our analysis is the following 
prediction. Under the conditions that particle number 
fluctuations are frozen after chemical freezeout, the fluc- 
tuations of non-conserved quantities, such as, e.g., mean 
Pt, are also preserved over time scales longer than colli- 
sional relaxation time r. The strength of this effect cru- 
cially depends on the ratio of the a screening mass mt 
at kinetic freezeout to A (see Eq. (19"T1) and Section IX D I) . 

In other words, while for the multiplicity fluctuations 
to be preserved after chemical freezeout the kinematic 
window of acceptance must be large enough [25l |26| , for 
the pr fluctuations to be preserved, additional condition 
is necessary: < A. We find that, e.g., at top SPS 
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energies, pt fluctuations can survive the hadronic rescat- 
tering at least half as well as the particle multiplicity 
fluctuations for < 300 MeV. 

The origin of this effect is the mixing between the crit- 
ical mode a and the conserved particle number density 
(see discussion after Eq. ([571) ). E.g., when A » m, the 
mode h\ , which is kept from relaxing by the particle 
number conservation, is almost the same as ct, Eq. ()62[) . 
The fluctuations of a involved in h\ must keep the mag- 
nitude they reached at chemical freezeout, contributing 
the term ~ l/m 2 into Eqs. (193)) and (|96|) . While multi- 
plicity fluctuations in Eqs. (|93|) are frozen, the pt fluc- 
tuations evolve, with contribution of the mode h\ a de- 
creasing with increasing m as the factor r m . 

At the same time, the fluctuations of a alone, with 
particle number fixed (i.e., obeying Eq. (1581) ). equilibrate 
on a short time scale, r , tracking the evolution of m. 
This equilibrated mode of fluctuations contributes 1 / m£ 
term into Eq. ((96]) . 

In this paper we focused on fluctuations of one particle 
species, treating the rest of the hadron gas as a heat bath. 
This simplification allowed us to follow the evolution of 
fluctuations analytically and expose the mechanism be- 
hind the "memory" effect in the most transparent way. 
This analysis could be generalized to the case of multiple 
particle species, carrying (different values of) the same 
conserved charge, as well as the case of multiple con- 
served quantities (baryon number, isospin, etc.). Taking 
into account fluctuations of conserved energy and mo- 
mentum would be necessary, for example, to obtain cor- 
rect vtl — y scaling of the diffusion coefficient [111 |43[ ■ 
We leave this to future work. 

We also neglected the effects of quantum statistics 
for simplicity. Although these are relatively small un- 
der realistic conditions (few percent, as estimated by 
mean occupation numbers Ai[fo\), this approximation 
could be removed. For the most part this would re- 
quire replacing the equilibrium distribution in Eq. (|10p 
with Bose-Einstein or Fermi-Dirac distribution and fac- 
tors fo in equations such as Eq. ([52} with /o(l ± fo)- 
This would also imply that the collision integral C [f] has 
Uehling-Uhlenbeck form [5(j. The influence functional 
method [HJ HH could be used to derive the correspond- 
ing equations. 

We would like to stress that, although we did use 
relaxation-time approximation to obtain more explicit 
formulas for non-zero modes A in Section IIX1 the re- 
sults pertaining to the "memory" effect, which rely on the 
properties of the zero mode Ao studied in Section IVIII1 
are valid beyond relaxation-time approximation. 

A numerical simulation of the stochastic equations (j!8[) 
should allow to take into account more detailed prop- 
erties of the heavy-ion collision evolution, such as in- 
homogeneity, anisotropy and flow. The evolution of 
the a mass m can be described self-consistently, using 
equation (|22p . conceptually reminiscent of nuclear mean- 
field approach [34| , or disoriented chiral condensate stud- 
ies [53L [54j. 



We also deliberately limited our analysis to linearized 
regime and focused on quadratic moments of fluctua- 
tions. The stronger singular behavior of higher moments 
of fluctuations makes them more attractive signatures of 
the QCD critical point [55[. A study of the higher-order 
moments would require generalization of the analysis to 
nonlinear equations such as ((5} and ([9j. 
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Appendix A: Notations 

C[f] - collision integral (01; 
D - diffusion constant ((72)) ; 
Do _ same, at g = 0, ((7T)) : 
E - "entropy matrix" ([23 ]) . ([23" ]> : 

/ - short for f(x,p;t), non-equilibrium distribution 
function ((4)); 

f a — equilibrium distribution function for given back- 
ground a ((TO)) : 

fo or /q - short for fo(p) or fo{p'), as above, for a = 0; 

G - operator ((27)) . acting on h gives "drift" terms in 
stochastic equations (fT8)). ((26)) ; 

g - coupling of particles to a ((14")) ; 

h - short for h(x,p;t) relative deviation of / from 

fa CCD) ; 

h - mean value of h ([56)) : 

h - deviation of h from h ([56)) : 

h - generalized vector (h, tt, a) ((24)) : 

h\ - eigenmode of G with eigenvalue A ([521 : 

h - convenient basis vector (|59l) : 

I - linearized collision integral (TTS)) ; 

JC - linear operator related to I by Eq. ([3T1) ; 

M or M(a) - particle mass for given a ([T]); 

M [h] - mean value of h ([49]) ; 

mo - vacuum mass of the field a ((TC)) . d 2 U(0)/da 2 ; 

to - thermal screening mass of a (| 1 7[) , ([22)) ; 

to - thermal pole mass (rest energy) of a quasiparti- 

cles ([79]), ([81)1: 
to c or TOk - screening masses of a at chemical or kinetic 

freezeout; 
N - total number of particles; 
no - equilibrium density of particles (|55|) : 
p or p' - particle momentum variable in /; 
Pt - the magnitude of the component of p transverse 

to beam axis; 
q - Fourier conjugate to x in h(x,p;t); 
Q - matrix of noise correlators l[2"5)) ; 
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r m - "memory" factor (|97|) ; 

- quadratic terms in the entropy ([2"Tjl ; 
T - temperature of the external bath ([20)) ; 
U - U(a) potential for a, ([1]); 

V - J d 3 x, 3-volume; 

V - evolution operator (j43|) : 
v - short for v(p), particle velocity (j5]); 
To - relaxation rate of a due to interaction with the 

external thermal bath only ()18bl) ; 
T - the full relaxation rate of a ([75]). (IM)) ; 
7 or 7' - 7(p) or 7(p')> relativistic factor ©; 
A 2 - Eq. {63) and also m 2 - m 2 {JJ]); 

V - (2^) 3 <5 3 (p-p'), ©; 

77 - 77(2;), Langevin noise in eq. (|18bj) : 
0\ - dual vector to h\ 
- see Eq. (JUD ; 



Ao - smallest eigenvalue of (|52|) ; 

i'p - momentum space distribution, (|37p ; 

£ - t;(x,p;t), noise in Boltzmann eqn. (|18a[) ; 

1 - noise vector 77, 0) (gU), O; 

7r - canonical momentum for a, (|18c[) ; 

X - matrix of correlators (h <S> h), (|4"Tj) ; 

S - deviation of S from equilibrium (|45p ; 

a - <j{x), scalar field, critical mode; 

<r - basis vector (|BU|) ; 

r - collisional relaxation time (174)) ; 

r e - fireball evolution time scale, Sec. IX Al 

%l) - ip(v 2 ), solution to Eq. (|55|) ; 

/ p - see Eq. ([8J; 

• - the scalar product is defined in Eq. (|2"5|) . 
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